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Abstract 

Gaussian sum-rules, which are related to a two-parameter Gaussian-weighted integral of a hadronic spectral 
function, are able to examine the possibility that more than one resonance makes a significant contribution to 
the spectral function. The Gaussian sum-rules, including instanton effects, for scalar gluonic and non-strange 
scalar quark currents clearly indicate a distribution of the resonance strength in their respective spectral func- 
tions. Furthermore, analysis of a two narrow resonance model leads to excellent agreement between theory and 
' phenomenology in both channels. The scalar quark and gluonic sum-rules are remarkably consistent in their 

, prediction of masses of approximately 1 GeV and 1.4 GeV within this model. Such a similarity would be expected 

' from hadronic states which are mixtures of gluonium and quark mesons. 
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1 Introduction 

Oh! 

^ . The multitude of scalar hadronic states with masses above 1 GeV PP is frequently noted as evidence that a qq nonet 
'T^ I is insufficient to accommodate these states, as would be expected from the existence of a (scalar) gluonium state. 
^ ■ Gaussian QCD sum-rules have been shown to be sensitive to the hadronic spectral function over a broad energy 
I range, and analysis techniques have been developed to exploit this dependence to determine how resonance strength 
. is distributed in the spectral function |21I31S- Thus Gaussian sum-rules provide a valuable technique for determining 
I how the quark and gluonium content is distributed amongst scalar hadronic states. 
The simplest Gaussian sum-rule (GSR) has the form j5J 



Go (s,t) = -L 



OO 




/exp 




to 





-p{t)dt , r>0 (1) 

TT 

and relates the QCD prediction Gq (J, t) to an integral of its associated hadronic spectral function p(t). The smearing 
of the spectral function by the Gaussian kernel peaked at t = s through the (approximate) region s — < t < 
s + provides a clear conceptual implementation of quark-hadron duality. The width of this duality interval is 
constrained by QCD since renormalization-group improvement of the QCD (left-hand) side of Q results in identifying 
the renormalization scale ly through ly^ = -y/r [2 E| ; therefore it is not possible to achieve the formal t ^ limit 
where complete knowledge of the spectral function could be obtained through 

lim Go (s, r) = -p{s) , s > to . (2) 

T— >0 TT 

The variable s in (Q, on the other hand, is unconstrained by QCD, and so the s dependence of Go (s, t) can be used 
to probe the behaviour of the smeared spectral function, and hence the essential features of p{t). 

An interesting feature of the GSR is its ability to study excited and ground states with similar sensitivity. For 
example, as s passes through t values corresponding to resonance peaks, the Gaussian kernel reaches its maximum 
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value. Thus any features of the spectral function strong enough to be isolated from the continuum will be revealed 
through the GSR. In this regard, GSRs should be contrasted with Laplace sum-rules 

oo 

i?(A2)=i Jexp(^~±y{t)dt , (3) 

to 

which exponentially suppress excited states in comparison to the ground state. 

In what follows, the original formulation of GSRs will be reviewed along with corresponding GSR analysis 
techniques |2l |3j . A GSR moments analysis for scalar gluonic and quark correlation functions will be presented in 
order to demonstrate that the associated spectral functions have a distribution of resonance strength. In addition, 
mass and relative coupling predictions specific to a double narrow resonance model will be determined [3 IHj ■ The 
major emphasis of this paper is to highlight the remarkable consistency between GSR mass predictions extracted 
from the scalar gluonic channel and from the scalar-isoscalar quark channel, a result indicative of the existence of 
hadronic states which are mixtures of quark mesons and gluonium. 



2 Formulation and Analysis Techniques for Gaussian Sum-Rules 

Gaussian sum-rules are based on QCD correlation functions of renormalization-group invariant composite operators 
Jix) 



n (g2) ^ijd^x e'«-- (0\T [J(x) J(0)] \0) , 



(4) 



which, in turn, satisfy dispersion relations appropriate to the asymptotic form of the correlator in question. For 
example, the scalar gluonic correlation function 11 (Q^) (see (|24ll and (jUJ) satisfies the following dispersion relation 
with three subtraction constants 



n (Q V n(0) - WO) - lQ^n"(o) ^ 4 ^ ^3 ^ 



Pit) 



dt 



(5) 



In general, the quantity p{t) is the spectral function appropriate to the quantum numbers of the given current, 
and it should be noted that in certain situations (such as the scalar gluonic correlation function) the subtraction 
constant 11(0) is determined by a low-energy theorem [S] (see (|29|) 'l. Undetermined dispersion-relation constants and 
field-theoretical divergences are eliminated through the GSR j3)^ 



Gfc(s,T) 



/ (s iA)'^-n(-.s - iA) - (s - iA)'=n(-s + iA) 



(6) 



where /c = — 1,0,1,... and where the Borel transform B is defined by 

,2\N / -1 \ N 



B = lim 



iv,A--.oo r(iv) 



dA2 



(7) 



The dispersion relation (jSJ in conjunction with definition © together yield the following family of GSRs 



Gfc(s,T) + 4- 

where the identity 
1 



'Attt 



exp _ n(G) = 



4t 



/47rr 



t'' exp 



to 



At 



-p{t) dt 

TT 



B 



A2 -f a 



4t 



At 



— exp — , n = 0,l,2,. 



(8) 



(9) 



has been used to simplify the phenomenological (right-hand) side of (jSJ as well as the term on the QCD side 
proportional to 11(0). Evidently the k — —1 sum-rule can only be defined in cases where there exists an appropriate 

^ This definition is a natural generalization of that given in |5] . To recover the original GSR, we simply let fe = in j^J . 
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low-energy theorem. Calculation of G'fe(s,r) is achieved through an identity relating Q to the inverse Laplace 
transform |S] 



B[fiA')] = -C-'[fiA')] 
where, in our notation, 



a+ioo 



(10) 



(11) 



with a chosen such that all singularities of / lie to the left of a in the complex A^-plane. Through a simple change 
of variables, the calculation of the GSR reduces to |2| 



Gfe(s,T) 



1 1 



/Attt 27ri 7ri+r2 



{-w)'' exp 



4r 



Il{w) dw 



which, by suitably deforming the complex contours Fi and can be recast in the form 



Gfc(s,T) 



1 1 



^Attt 27ri Jr^+r 



{-w)'' exp 



- (s + w) 



21 



4t 



Il{w) dw 



(12) 



(13) 



where the trajectories Fi, r^, and are all depicted in Figure^ Further simplification of (|13|) requires a specific 
correlator 11 i.e. or (|77|l . 




Figure 1: Contours of integration Fi + r2 and Fc + F^ defining the Gaussian sum-rule (see 11211 and I13II respectively). The 
wavy line on the negative real axis denotes the branch cut of Il{w). 



Next, we impose a fairly general resonance(s) plus continuum model 

p{t) = phad(^) + e{t^so) imnQCD(t) 



(14) 



where sq represents the onset of the QCD continuum. The continuum contribution of H14|l to the right-hand side 
of © is 



Gr'(s,r,so) 



/47rT 



t'^ exp 



4r 



-imnQ™(0 dt , 



and is combined with Gk (s, t) to obtain the total QCD contribution 
gQ^° {s, t, so) EE Gfe is, t) - is, T, so) , 



(15) 



(16) 
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resulting in the final relation between the QCD and hadronic sides of the GSRs 

S2\ 1 /-oo ^-{s-ty 



is, r, .0) + 4 exp (^) n(0) = ^ 



t'' exp 



4t 



-p'^^^(t) dt 

TT 



Original studies involving the GSRs exploited the diffusion equation 



dr 



(17) 



(18) 



which follows from © for fc = 0. In particular, when p^^^{t) (see (I14II ') is evolved through the diffusion equation p8|l . 
it only reproduces the QCD prediction at large energies (r large) if the resonance and continuum contributions are 
balanced through the fc = member of the finite-energy sum-rule family |5j 



1 



So 



TT 

to 



(19) 



An additional connection between the finite-energy sum-rules ()19|l and the GSRs can be found by integrating 
both sides of (|17|) with respect to s to obtain 



oc oo 

y" G^^°(s,r,so) dJ + 4-in(0) = i j t''p'^-''{t)dt , 



(20) 



to 



indicating that the finite-energy sum-rules (|19|l are related to the normalization of the GSRs. Thus the information 
independent of the finite-energy sum-rule constraint following from the diffusion equation analysis |S| is contained in 
the normalized Gaussian sum-rules (NGSRs) |2] 



(21) 
(22) 



^ Gg^"(g,r,.o)+4-r^exp(^)n(0) 
' MQo^°(r,.o) + 4-in(0) 

oo 

^fe,n(T, So) = y s"Gfe(s,r, So) ds , n = 0,1,2,... , 



which are related to the hadronic spectral function via 



7VQ^°(s,r,so) = 



1 rOO ,1. 

J= . rexp 



p'^^^(t) dt 



(23) 



3 Gaussian Sum- Rules for Scalar Gluonic and Non-Strange Scalar Quark 
Currents 

At leading order in the quark mass for nf flavours, scalar hadronic states can either be probed through the correlation 
function of the (renormalization-group invariant) gluonic current 



n,(g2) ^i / d'xe"'-^{0\T[J,{x)J,m\0) , 



Ja{^) = --^P{a)G%{x)GlA^) 



(24) 
(25) 



where 



11 1 51 19 



(26) 
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or through the correlation function of / = 0, 1 (non-strange) quark currents 



U,{Q^)^i J d^xe'^--(0\T[Jg{x)J,{0)]\0) , 
Jg{x) = mq \u{x)u{x) + {-lY2{x)d{x)\ /2 . 



(27) 
(28) 



The quark mass factor niq — {niu + md) /2 in H28|l results in a renormalization-group invariant current. Although 
both the gluonic and (/ — 0) quark correlation functions are probes of scalar hadronic states, those which have a 
more significant overlap with the gluonic current will predominate in (|24|l . while those states which are dominantly 
of a quark nature are more significant in 1)2 7|l . A mixed state with substantial gluonic and quark components {i.e. a 
state that overlaps with both the gluonic and quark currents) should self-consistently appear in an analysis of both 
correlation functions. In particular, prediction of mass-degenerate states from the QCD sum-rule analysis associated 
with these two correlation functions is evidence for interpretation as a mixed state. 
In the scalar gluonic channel, the low-energy theorem (LET) 



ng(0) = lim ..g 



ng(Q' 



Po 



(29) 



allows construction of the fc = — 1 GSR. The significance of instanton contributions in the overall consistency of 
the LET-sensitive fc = — 1 sum-rule and the LET-insensitive fc > sum-rules was first demonstrated for Laplace 
sum-rules [3IH1- A similar consistency is observed for the Gaussian sum-rules, but theoretical uncertainties are better 
controlled in the fc > GSR 0, and hence this paper will focus on the fc = GSRs for both the quark and gluonic 
channels. 

The QCD correlation functions (|24|l and (|27|l contain perturbative, condensate, and instanton contributions. For 
the scalar gluonic case, substituting H24|l into H13|l and performing the relevant integrals yields the QCD prediction 

Gq^^ (s, r, So) for the fc = GSR to leading order in the quark mass 



Gi'\s,T, so)=- 




4t 



(ao - TT 02) + 2ailog ( — 



3a2 log ( ^ 



dt 



1 




exp 





So 



At 



: exp 



— s 
17 



CO {Oe 



(30) 



exp 



At 



The perturbative coefficients in (|3()|l are given by 

247.480 



ao 



1 



659 a 



ai 









65.781- 




4 TT 



0,2 



-10.1250 f- 
Vtt 



(31) 



as obtained from the three-loop MS calculation of the correlation function in the chiral limit of = 3 flavours "5". 
The condensate contributions in H30() involve next-to-leading order (1011 contributions^ from the dimension four gluon 
condensate {Jg) and leading order ,11) contributions from gluonic condensates of dimension six and eight 



(Ce) — (gfabcG'^^GlpGp^) 

(Os) = lA({aU,G^ppGlp)') - UaUcGl^G'p^) 



bn = An— 



Co = 87r- 



175 a 



b,=-97T - 



TT 



^The calculation of next-to- leading contributions in IIUI have been extended non-trivially to n^^ 
basis has been changed from (aG^) to (Jg). 



(32) 
(33) 

(34) 

0, and the operator 
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The remaining term in the GSR (|30|l represents instanton contributions obtained from single instanton and anti- 
instanton |12| (i.e. assuming that multi- instanton effects are negHgible |13|) contributions to the scalar gluonic 
correlator [7ll51 lll[ll4| within the liquid instanton model |15| parameterized by the instanton size p and the instanton 
density ric. The quantities J2 and Y2 are Bessel functions in the notation of 

As a result of renormalization-group scaling of the GSRs |21 E] , the coupling in the perturbative and condensate 
coefficients f (|S1T) and is implicitly the running coupling at the scale v"^ = yjr for n/ = 3 in the MS scheme 

a{v'^) _ 1 /3ilogL 1 



/3oi (/3oL)' (/3oi)' 



^ = log ( ^ 



(log2L-logL-l)+/32] 

3863 



(35) 



384 



with AjY^ sa 300 MeV for three active flavours, consistent with current estimates of q;(Mt-) 

Analogous to the scalar gluonic case, substitution of \27\i into H13|l provides the QCD prediction of Gq'''' (s, t, sq) 
for the 7 = 0,1 scalar quark currents. To leading order in the quark mass, we find P] 



Gi'^ {S,T, so) 



1 



3ml 



exp 



exp 



-jt-s) 

At 



M 1 + ?- ) -2-<log( ^ 



3 TT 



4t 



1 



-(-1) 



2^ 
t exp 



4ry'7rr 



4r 



Ji U%/7j Yi fpVij di 



dt 



(36) 



Again, renormalization-group improvement implies that both and a are implicitly running quantities at the scale 
v"^ = ^/t as given by and the (two-loop, n/ = 3, MS) expression 



290 1 256 log L 
729 L ^ 729 L 



L = log 



(37) 



where rhq is the renormalization-group invariant quark mass parameter. The perturbative contributions in H36II are 
the 71/ = 3 two-loop results obtained from |17j , and the instanton expressions are obtained from |18j . The condensate 
contributions are leading-order results obtained from ^H], and are defined by the quantities 



{ClOD = \ (m.qq) + ^ (a^G^) 



(38) 



and 




u,d,s I 



(39) 



The vacuum saturation hypothesis ^Hl in the SU (2) limit {uu) = (dd) = (qq) provides a reference value for (0| 



(QO^) = -/„,-a, {(qqf) = -/,,1.8 x lO-^GeV^ 



(40) 



where the quantity fyg parameterizes deviations from vacuum saturation where fvs ~ 1. 

The GSRs (|30|) and (|36|l exhibit some interesting qualitative features. For example, the condensate contributions 
decay exponentially with the Gaussian peak-position s, emphasizing that these contributions have a low-energy 
origin. Also, the explicit factor of / appearing in the instanton contributions in the quark scalar channel are a 
non-perturbative source of isospin symmetry-breaking. 

Before proceeding with an analysis of the GSRs, the QCD input parameters must be specified. We assume that 
(Jg) « (aG^) and then employ the (central) value from 



{aC^^C^"") = (aG') = (0.07 ± 0.01) GeV^ 



(41) 
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The dimension six gluon condensate H32|) can be related to this value of {aG'^) using instanton techniques (see |ll[ll8p 



(Oe) = (0.27GeV^)(aG2) . (42) 

Further, by invoking vacuum saturation in conjunction with the heavy quark expansion, the authors of {2(Jj have also 
related the dimension eight gluon condensate to (aG^) through 

(08) = ^((aG2))' . (43) 
In addition, the dilute instanton liquid (DIL) model |15| parameters 

71c = 8.0 X 10"'' GeV^ , p= — GcV-^ (44) 
' 0.6 ^ ' 

will be employed. Finally, we use fvs = 1.5 as a central value to accommodate the observed deviations from vacuum 
saturation in the dimension six quark condensates H40|l |21) . Note that knowledge of the quark mass parameter is 
not needed for an analysis based on the normalized GSRs since ruq provides a common prefactor in (|36|l . 



Analysis of the Gaussian Sum-Rules for Scalar Gluonic and Quark 
Currents 



In the single narrow resonance model 

the /c = NGSR ^ becomes 

1 \-{s^m^f 



7VQ^° (s,r,so) 



: exp 



4r 



(45) 



(46) 



The prediction of the mass m is obtained by optimizing the parameter sq so that the left-hand side of l|46|) has a 
maximum as a function of s (s peak position) independent of r as required by the properties of the right-hand side 
of H46|l |2j. The r-stable s peak for this optimized sq then provides the prediction of the resonance mass m. 

The integrity of this procedure has been demonstrated for the vector-isovector currents which probe the p meson, 
resulting in a predicted p mass of 750 MeV and superb agreement between the phenomcnological and QCD sides 
of the NGSR However, the same procedure applied to the scalar gluonic and scalar quark channels leads to 
significant disagreement between the two sides of the NGSR as illustrated in Figure [2 ^ J .'^ In both the scalar 
gluonic and scalar quark sum-rules, the single resonance model is larger than the theoretical contribution at the 
peak and underestimates the theoretical contribution in the tails. Since both the theoretical and phenomcnological 
contributions are normalized, this implies that the QCD prediction is broader than the phenomenological model. 

The following second-order moment combination 



M, 



0,2 



(47) 



provides a quantitative measure of the width of the GSRs. In the single narrow resonance model we find cr^ — 2t, 
and hence a significant deviation of the QCD prediction from this result indicates a failure of the phenomenological 
model to adequately describe a particular channel's hadronic content. Figure O illustrates that > 2t, indicating 
that a phenomenological model with distributed resonance strength is necessary j3 13] . 

A phenomenological model with two narrow resonances of mass mi and m2 results in the NGSR 



^^0 (s,T, So) 



ri 



: exp 



— s — m. 



4t 



: exp 



4r 



(48) 



^The r-stability analysis of the s peak is examined in the range 2GeV* < t < 4GeV* where the perturbative series is reasonably 
convergent while maintaining a Gaussian resolution of typical hadronic scales. 
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Figure 2: Comparison of the theoretical prediction for the normahzed GSR N^f' (s,r, Sq) with the single narrow 
resonance phenomenological model for the optimized value of the continuum sq. The r values used for the three 
pairs of curves, from top to bottom in the figure, are respectively r = 2.0 GeV, r — 3.0 GeV'', and r = 4.0 GeV*. 
A qualitatively similar agreement between the single narrow resonance model and the QCD prediction exists for the 
/ = 0, 1 scalar quark channels. 



where ri + 7'2 = 1 describes the relative strength of the two resonances contributing to the spectral function. In 
terms of the parameters {r^y, z} defined by 

r = ri — r2 , y = ml — m\ , z — m\ + , (49) 

the second-order moment combination (|47|l resulting from the left-hand side of ()48|l 

CT^ - 2r = (^i _ 7-2) > (50) 

naturally results in a broader distribution than the single resonance model as suggested by the QCD result. 

Analysis of the double narrow resonance model is substantially more complicated than the single narrow resonance 
case. In particular, the s peak position (denoted by Speak) develops t dependence which is well-described by [5] 

B C 

Speak (r. So) = ^ H 1- ^ (51) 

and hence so is optimized to obtain the best description of this r dependence. After optimization of sqj the parameters 
in the phenomenological model can be determined from the moment combinations 

. = 2^ + ^^ (52) 
Mo cr2 - 2r ^ ' 



(54) 



^Ai + 4(a2 - 2r)3 

where the third-order moment combination representing the asymmetry of the distribution, is defined by 

A2 = ^-3^^ + 2r^V . (55) 

Mo,o Mo.oMo,o \Mo,qJ 

This procedure for optimizing sq and determining the resonance parameters has been confirmed by a more numerically- 
intensive multi-parameter fit of sq and the resonance parameters 22^ . 

The resonance parameters resulting from this analysis shown in Table ^[2121 illustrate a remarkably consistent 
scenario of a 1 GeV and a 1.4 GeV state coupled to both the gluonic and non-strange 7 = quark currents, with 
the heavier state more strongly coupled to the gluonic operators. This consistency of the mass predictions in the 
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T (GeV*) 



Figure 3: Plot of for the theoretical prediction (dotted curve) for the scalar gluonic GSR compared with cr^ = 2r 
for the single-resonance model (solid curve) using the optimized value of the continuum. A qualitatively similar 
result exists for the / = 0, 1 scalar quark channels. 



Sum-Rule 


mi (GeV) 


ma (GeV) 


ri 


7-2 


so (GeV^) 


gluonic 


0.98 


1.4 


0.28 


0.72 


2.3 


quark / = 


0.97 


1.4 


0.63 


0.37 


2.6 


quark / = 1 


1.4 


1.8 


0.57 


0.43 


3.9 



Table 1: Analysis results from scalar quark and gluonic Gaussian sum-rules in the double narrow resonance model. 

two channels is precisely what is expected for hadronic states which are mixtures of gluonium and quark mesons. 
The results in the 1=1 scalar quark channel support the interpretation of the ao(1450) as the lightest state with a 
dominant coupling to the scalar quark currents. 

The double narrow resonance model results in excellent agreement with QCD as illustrated in Figure 01 More 
complicated resonance models which extend the double narrow resonance model to introduce resonance widths do 
not improve the quantitative agreement with QCD exhibited in the figure. 

In summary, analysis of the GSRs for the scalar gluonic and / = non-strange quark scalar currents exhibit 
a remarkable similarity in their mass predictions within a double narrow resonance model, a result indicative of 
the existence of hadronic states which are mixtures of quark mesons and gluonium. Although the effect of QCD 
uncertainties arising from a 15% variation in the DIL parameters (|44f) and the d = A gluon condensate (|41|l . as well 
as variations of the vacuum saturation parameter within the range 1 < /„s < 2 correspond to an uncertainty in the 
Tabled mass parameters of approximately 0.2 GeV, we note that the mass splitting of 0.4 GeV between the states 
is remarkably stable [H]- This QCD evidence for mixed states with a mass splitting of 0.4 GeV provides valuable 
information for interpretation of the known /q resonances. 
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